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Chapter  1  Introduction 

As  the  speed  of  transistors  has  risen  in  recent  years,  acquiring  accurate  models  for  the 
passive  devices  used  in  integrated  circuits  has  become  a  greater  challenge.  Current  device  speeds 
are  in  the  tens  of  GHz,  and  at  these  frequencies  the  parasitic  effects  of  inductors  can  become 
difficult  to  predict.  At  GHz  level  frequencies,  monolithic  inductors  suffer  from  severe  parasitic 
effects.  The  relatively  large  size  of  the  inductor  (hundreds  of  microns),  and  the  high  dielectric 
constant  of  the  substrate,  cause  the  inductor  to  experience  significant  capacitive  coupling.  The 
skin  effect,  along  with  the  fields  coupling  into  a  lossy  substrate,  causes  the  structures  to  dissipate 
a  significant  amount  of  power.  Additionally,  the  effect  of  eddy  currents  must  be  understood  as 
they  too  result  in  power  loss.  These  parasitic  effects  result  in  a  limited  peak  Q.  Since  many  RF 
circuits  such  as  LNA’s  and  VCO’s  require  inductors  with  a  high  Q,  it  is  critical  for  the  circuit 
designer  to  accurately  model  these  parasitics. 

At  present,  the  best  way  to  model  the  parasitic  effects  of  a  monolithic  inductor  is  to  apply 
data  fitting  techniques  to  a  lumped  circuit  model,  or  to  use  a  numerical  EM  solver.  Data  fitting 
techniques  require  the  inductor  to  be  fabricated  multiple  times  which  is  both  time  consuming  and 
expensive.  EM  solvers  such  as  FEM,  MoM,  and  FDTD  can  be  both  time  consuming  and 
numerically  unstable  when  simulating  electrically  small  structures  (inductors  are  much  smaller 
than  X/20).  Although  new  schemes  have  been  developed  to  handle  electrically  small  structures, 
numerical  techniques  are  still  time  consuming. 

There  are  techniques  in  the  literature,  such  as  in  [17]  that  use  microstrip  transmission  line 
approximations  to  develop  general  expressions  for  the  inductance  and  Q  of  spiral  inductors. 
Other  techniques,  like  that  of  [1],  [6],  [8],  [15],  and  [16],  involve  the  development  of 
approximated  empirical  or  semi-analytical  formulas.  Unfortunately,  as  mentioned  by  [11],  these 
formulas  can  sometimes  give  results  with  an  error  on  the  order  of  20%. 

In  this  report,  we  begin  in  Chapter  2  with  a  summary  of  the  Partial  Element  Equivalent 
Circuit  (PEEC)  distributed  circuit  model  for  a  spiral  inductor.  Chapter  3  will  deal  with  the 
derivation  of  the  frequency  dependent  resistance  model  that  accounts  for  skin  effect  and  eddy 
currents.  In  Chapter  4,  we  review  the  approximations  made  to  simplify  the  Partial  Element 
Equivalent  Circuit  (PEEC)  model  into  a  simple  lumped  pi-model.  In  Chapter  5,  we  will  see  that 
the  lumped  model  can  be  evaluated  fairly  quickly  and,  as  will  be  demonstrated,  with  great 
accuracy.  We  will  summarize  our  results  in  Chapter  6. 
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Chapter  2  Summary  of  the  PEEC  Distributed  Circuit 

The  PEEC  technique  has  been  thoroughly  developed  by  [4]  as  a  tool  by  which  Maxwell’s 
equations  are  used  to  construct  a  distributed  circuit  to  model  interconnects  in  digital  and 
microwave  circuits.  As  the  construction  of  the  PEEC  distributed  circuit  for  inductors  is  discussed 
in  great  detail  in  [1],  [4]  -  [6],  and  [9],  it  is  only  briefly  summarized  here,  with  reference  to  [1 1]. 
A  flowchart  showing  the  progression  from  Maxwell’s  equations  (with  a  quasistatic 
approximation)  to  a  distributed  model  is  given  in  Figure  1. 


Figure  1  -  Flowchart  showing  the  Distributed  Model  Derivation 
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2.1  Review  of  EM  Equations 

Using  Faraday’s  law  in  point  form  for  the  time  periodic  case  concurrently  with  the 
constitutive  relations  for  linear  isotropic  materials,  one  can  arrive  at  the  following  expression  for 
the  electric  field. 


E  =  -jcoA-V  (f)  (2.1) 

Taking  the  divergence  of  both  sides  of  equation  (2.1)  and  applying  the  Lorenz  gauge  to  the  result 
will  give 

V20  +  CO2jUS<fi  =  —^-  (2.2) 

£ 

where  (j)  is  the  scalar  electric  potential. 

In  a  similar  manner  beginning  from  Ampere’s  law,  we  arrive  at  the  magnetic  equivalent 
to  equation  (2.2) 

V2A  +  6)2jusA  =  -juJ,  (2.3) 

where  A  is  the  vector  magnetic  potential.  Since  the  devices  we  are  studying  are  electrically  very 
small,  we  can  apply  a  quasi-static  (co  ~  0)  assumption  to  equation  (2.2)  and  equation  (2.3)  which 
gives 

VV  =  -—  (2.4) 

£ 

V2A  =  -juJ  (2.5) 

Equation  (2.5)  is  based  upon  the  assumption  that  the  substrate  has  very  low  loss.  If  we  want  to 


take  into  account  losses  in  the  substrate,  we  can  use ss  =  £(]£rs  -  j—1-  in  equation  (2.3)  to  get 

co 


(within  the  substrate) 


V  A  +  co  /Jq 


.  cr 

S0£rs  ~~J  — 

v  CO  j 


AU  V2A-  jcoo-sju0A  =  ~/uJ , 


(2.6) 


where  the  co2  term  is  ignored  because  we  are  using  a  quasi-static  approximation.  The  solution  to 
equation  (2.5)  and  equation  (2.6),  with  the  appropriate  boundary  conditions,  is  expressed  in 
tenns  of  the  magnetostatic  ( Gm )  Green’s  function. 

A  =  f  G  Jdv  (2.7) 

J  v 
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The  solution  for  Gm  for  the  lossless  case  (equation  (2.5)  )  is  simple  because  it  does  not  involve 
the  permittivity.  Since  the  permeability  of  each  layer  of  our  substrate  is  the  same  as  free  space, 
Gm  is  simply  the  free  space  Green’s  function  without  the  e~jBR  term,  namely, 

G"W  =  tV  (2-8) 

AnR 

where  R  is  the  distance  between  the  observation  and  source  points.  The  solution  to  equation  (2.6) 
(lossy  case)  is  more  involved  and  requires  the  application  of  appropriate  boundary  conditions  to 
solve  for  a  multi-layer  Gm.  That  case  will  not  be  considered  in  this  report.  The  solution  to 
equation  (2.4),  with  the  appropriate  boundary  conditions,  can  be  expressed  in  terms  of  the 
electrostatic  ( Ge )  Green’s  function, 

(f>  =  f  G  pdv .  (2.9) 

J  V 

As  mentioned  previously,  the  inductors  we  are  studying  are  electrically  very  small, 
making  the  use  of  static  Green’s  functions  appropriate.  Following  the  method  outlined  in  [5]  and 
[10],  Ge  for  the  multi-layer  substrate  shown  in  Figure  2  is  solved  directly.  The  details  of  the 
derivation  of  the  multilayer  Green’s  function  Ge(r,r')  are  given  in  Appendix  A. 

h»E  =  0 

"A 

METAL  3 
METAL  2 
METAL  1 
OXIDE 

SILICON  BULK 
GROUND  PLANE 

/ 

hxE  =  0 


m 

n 


Figure  2  -  Profile  of  multi-layer  substrate 

As  will  be  discussed,  equation  (2.7)  and  equation  (2.9)  fonn  the  foundation  for  the  system  of 
equations  that  will  be  used  to  construct  the  PEEC  distributed  circuit. 
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2.2  Discretizing  the  Inductor 


In  general,  the  current  and  charge  on  a  structure  can  be  discretized  any  arbitrary  way.  In 
other  words,  the  charge  and  current  can  be  expressed  in  terms  of  a  set  of  basis  functions  that  is 
most  appropriate  for  the  geometry.  However,  for  simplicity,  as  done  in  [1 1  ]-[  1 3],  the  monolithic 
inductors  are  discretized  into  a  set  of  series  rectangles  as  illustrated  in  Figure  3.  The  coordinate 
system  we  use  is  shown  in  Figure  4. 


Figure  3  -  Discretization  for  1.5  turn  inductor 


iZ 


Free  Space 

Z=h3 
Z=h2 
Z=h-, 
Z=0 


Substrate  Layer  1 


Substrate  Layer  2 


Substrate  Layer  3 


b) 


Figure  4  -  Coordinate  system  used:  (a)  top  view  and  (b)  side  view  of  substrate 
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The  current  and  charge  is  then  decomposed  into  a  summation  of  pulse  basis  functions  as 
done  in  [1 1]-[13].  The  volume  current  density  is  given  as 


J,  = 


Iib{x,y,z)iJj 


(2.10) 


where  b(x,y,z)i  is  a  pulse  function  with  a  magnitude  of  1  when  the  point  (x,y,z)  is  inside  rectangle 
i  and  is  equal  to  zero;  otherwise,  .7  is  a  unit  vector  that  gives  the  direction  of  current  flow  and  is 


always  perpendicular  to  the  cross-section  into  which  it  flows,  whose  area  is  denoted  by  at.  /,■  is  a 
constant  and  represents  the  magnitude  of  the  current.  Similarly,  the  volume  charge  density  in 
each  segment  is  given  as 


Qib(x,y,z)i 


(2.11) 


where  b(x,y,z)i  is  the  same  as  in  equation  (2.10),  v,  is  the  volume  of  rectangle  i,  and  0,  is  a 
constant  that  represents  the  total  charge  in  the  volume  v,.  The  total  current  and  charge  in  the 
inductor  is  a  summation  of  the  individual  basis  functions.  This  gives  the  total  current  and  charge 
density  as 


N 

p=Tjp< 

i= 1 


(2.12) 

(2.13) 


2.3  Formulation  of  the  PEEC  System  of  Equations 

It  is  now  possible  to  construct  the  PEEC  distributed  circuit  with  the  equations  we  have 
derived  up  to  this  point.  Plugging  equation  (2.7)  into  equation  (2.1)  and  putting  all  tenns  to  the 
right  side  will  give 

0  =  E  +  ja>  [  G  Jdv  +  V0 .  (2. 14) 

J  v 

We  can  replace  J  in  equation  (2.14)  with  equation  (2.10)  and  equation  (2.12)  to  arrive  at 

N  I 

o  =  E  +  jcoYj  —  £  Gmb(x,  y,  z).  J^v  +  Wfi.  (2.15) 

j= i  aj 
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By  defining  the  testing  function  to  be  the  same  as  the  basis  function  from  equation  (2.10),  we  can 
apply  Galerkin’s  method  and  take  the  inner  product  of  both  sides  of  equation  (2.15)  with  the 
testing  function  to  arrive  at 


0  =  1  +  |  J  G„{jpjl)dvJdvJ+\ 

JV:  n  H  V;  JV;  \  J  J  ■'  ■' 


»  I 


dv,  . 


(2.16) 


7=1  apj 


From  the  boundary  condition  on  a  good  conductor,  it  can  be  shown  that  the  first  tenn  on 
the  right  side  of  equation  (2.16)  represents  the  resistance  ( l/o )  of  the  inductor. 


-  -  J, 

E  •Jj  =—. 
a 


(2.17) 


We  call  the  double  integral  in  the  second  term  of  equation  (2.16)  My,  the  partial  inductance 
matrix,  namely, 

=  —  £  £  G,n  (  j,  •  Jj )  dvjdv, ,  (2- 1 8) 


Mi; 


aflj 


where  Gm  was  defined  in  equation  (2.8).  Plugging  equation  (2.8)  into  equation  (2.18),  we  get 


M,  =-£>- 

4  naa 


r  r  r  r  dl  ■  •  l//- 

-  — - - rdaJa,. 

J  Cl;  J/,  Jdj  Jlj  _  J4  \  d 


l  J 


(2.19) 


The  analytical  solution  to  equation  (2.19)  can  be  very  complicated.  The  techniques  outlined  in 
[2],  [6],  [7],  and  summarized  in  Appendix  C,  offer  approximate  solutions  to  equation  (2.19)  that 
give  reasonable  accuracy. 

Finally,  a  finite  difference  approximation  can  be  applied  to  the  third  term  of  equation 
(2.16)  and  simplify  the  equation  to 

1  N 

°=-J  l^Jdv:  ■  \F.  (2.20) 

Ui  V'  7=1 

By  definition,  AVt  is  the  change  in  voltage  from  one  end  of  rectangle  i  to  the  other  end.  As  in 


[11]-[13],  we  define  AVj  as 


A  Vi=<j)(x%  y*  ,zl)~  0(x; ,  y~  ,z~), 


(2.21) 


where  0(x(+,y,+  ,z;+)  denotes  the  potential  at  the  boundary  between  rectangle  i  and  rectangle  i-1. 


Similarly,  r/)(x,  ,  yt  ,  z;.  )  denotes  the  potential  at  the  boundary  between  rectangle  i  and  rectangle 
i+1.  Expressing  equation  (2.20)  in  matrix  form  gives 
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[A  V]  =  ([R]  +  ja)[M])[l].  (2.22) 

Whereas  the  system  of  equations  expressed  in  equation  (2.22)  is  based  on  equation  (2.1) 
and  equation  (2.7),  a  second  system  of  equations  is  generated  from  equation  (2.9)  and  equation 
(2.23),  the  continuity  equation. 

V  •/  =-jcop .  (2.23) 

Similar  to  the  application  of  Galerkin’s  method  used  in  arriving  at  equation  (2.16),  we  obtain 
equation  (2.24). 

Qj 
IL~ 

7=1 


V,=-\  t,dvl=±\  \ 

V  "  V,-  JV;  J  \ 


J  ViVj 


Gedvjdvr 


(2.24) 


Here,  Vt  is  defined  as  the  average  potential  on  rectangle  i.  If  we  assume  that  the  potential 
decreases  linearly  along  the  length  of  the  rectangle,  then  Vj  is  approximately  the  potential  at  the 
center  of  the  rectangle. 

In  using  this  technique  for  an  N  rectangle  structure,  there  are  two  variables  that  are  not 
assigned  to  rectangles,  namely  Vo  and  Vv+/-  However,  the  input  and  output  ports  of  the  inductor 
are  located  at  the  input  boundary  of  rectangle  1  and  the  output  boundary  of  rectangle  N 
respectively.  Therefore,  Vo  and  Vn+i  can  be  assigned  to  the  applied  port  voltage  sources. 

Qi  is  the  total  charge  present  on  rectangle  j  and  can  be  found  by  integrating  the  charge 
density  over  the  volume  of  rectangle  j. 

Qj=\  pdVj  (2.25) 

By  solving  equation  (2.23)  for  p  and  inserting  the  result  into  equation  (2.25),  substituting  the 
current  J  with  a  summation  of  pulse  basis  functions,  we  get 

(v  •I,Ji  )b(x,y,z). 


Qj=-\  r 

JCO  Jvj  t? 


-dvj. 


(2.26) 


The  divergence  is  taken  on  the  volume  current  density  in  rectangle  j.  The  result  is  a  derivative 
with  respect  to  the  direction  in  which  current  is  flowing.  Again,  a  finite  difference  approximation 
can  be  made  giving 


A/ 

Qj=  — 

JO) 


(2.27) 
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Here  AIj  is  the  difference  between  currents  at  the  input  and  output  boundaries  of  rectangle 
j.  Current  at  the  boundaries  is  assumed  to  be  the  average  of  the  two  rectangles  sharing  this 
boundary.  Here  we  let  70  =  Ix  and  IN+l  =  I N  .  By  substituting  (2.27)  into  (2.24),  we  obtain 


2.4  Creating  the  Distributed  Circuit  from  the  System  of  Equations 

Equations  (2.22)  and  (2.31)  are  in  the  form  of  Kirchoffs  voltage  law  (KVL)  and 
Kirchoff  s  current  law  (KCL),  respectively.  From  (2.31)  we  construct  the  circuit  shown  in  Figure 
5  using  KCL  with  some  algebraic  manipulation  as  in  [1 1  ]-[  1 3]  so  that  the  topology  of  the  circuit 
agrees  with  the  properties  of  the  capacitance  matrix  outlined  in  [14].  The  capacitance  of 
rectangle  i  to  ground  C„  can  be  written  as 

N 

Q=ZC/  (2.32) 

where  c,y  are  the  individual  terms  of  the  [C]  matrix  for  an  N  rectangle  system. 


Figure  5  -  Equivalent  KCL  for  three  series  rectangles 
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Also,  the  direct  capacitance  between  two  rectangles  i  and  j  can  be  written  as 

CV=~CV.  (2.33) 


This  leads  us  to  construct  the  circuit  so  that  the  current  to  ground  from  rectangle  i  is  written  as 
jcoCuVi  and  the  current  from  rectangle  i  to  j  is  written  as  jcoCj(Vr  Vj).  An  equivalent  circuit  of 
Figure  5  is  re-drawn  with  capacitors  in  Figure  6. 


Figure  6  -  Equivalent  KCL  with  capacitors 

To  complete  the  derivation  of  the  PEEC  distributed  circuit,  we  use  equation  (2.22)  to 
model  the  ohmic  loss  and  magnetic  coupling  within  the  metal  of  the  inductor.  Utilizing  the  KVL 
form  of  equation  (2.22),  the  circuit  in  Figure  7  is  developed. 


I 


Figure  7  -  KVL  circuit  for  one  rectangle 
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We  combine  the  circuit  of  Figure  7,  which  represents  the  magnetic  coupling  and  power 
loss  in  the  conductor,  with  the  circuit  of  Figure  6,  which  models  the  electrostatic  coupling 
between  different  conductive  segments  and  between  the  conductor  and  ground.  The  resulting 
distributed  model  is  formed  in  Figure  8. 


R3  M33 


Figure  8  -  PEEC  distributed  circuit  for  three  series  rectangles 

The  capacitors  in  Figure  4  are  complex.  This  is  related  to  the  fact  that  Ge  is  complex  and 
accounts  for  loss  in  the  substrate.  As  a  result,  each  partial  capacitance  is  represented  by  an  ideal 
capacitor  in  parallel  with  a  resistor.  The  ideal  capacitance  is  simply  the  real  part  of  the  partial 
capacitance.  The  parallel  resistance  is  simply  given  by 

R,J=  7  TV-  (2.34) 

c o[unag[C .JJ 

From  [5]  we  know  that  the  imaginary  part  of  Q  is  approximately  inversely  proportional  to 
frequency,  which  tells  us  that  Ry  is  approximately  constant  with  respect  to  frequency.  The  terms 
Rij  above  are  not  to  be  confused  with  the  terms  of  [A],  which  is  discussed  in  Chapter  3  and  deals 
with  the  frequency  dependent  resistance. 

The  application  of  the  distributed  model  to  the  inductor  is  simple.  It  is  observed  in  Figure 
3  that  the  rectangles  comprising  the  inductor  are  in  series  from  port  1  to  port  2,  allowing  us  to 
add  the  distributed  models  in  series.  Thus,  for  an  inductor  of  twelve  rectangles,  we  would  have 
four  of  the  circuits  from  Figure  8  in  series. 
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Chapter  3  Calculating  the  Resistance  Matrix 

In  Chapter  2,  the  current  distribution  is  assumed  to  be  uniform  inside  the  volume  of  the 
rectangle.  However,  the  current  will  be  non-uniform  due  to  the  skin  effect  and,  in  the  case  of  the 
spiral  inductor,  eddy  currents.  While  a  uniform  current  approximation  is  reasonably  adequate  for 
the  calculation  of  the  partial  inductance  and  partial  capacitance  matrices,  it  is  not  adequate  to 
calculate  the  resistance  matrix. 

3.1  Eddy  Currents  in  a  Spiral  Inductor 

We  begin  with  a  discussion  of  the  derivation  of  [/?].  The  resistance  matrix,  [/?],  is  a 
diagonal  matrix  whose  elements  represent  two  important  effects:  the  amount  of  power  dissipated 
in  each  rectangle  due  to  skin  effect  and  eddy  currents,  as  well  as  the  effects  due  to  the  internal 
inductance  of  each  rectangle.  To  calculate  the  total  power  dissipated  in  each  rectangle,  the  total 
current  in  each  rectangle  is  first  calculated.  As  an  approximation,  the  total  current  in  a  rectangle 
is  expressed  as  a  sum  of  an  excitation  current  and  an  eddy  current.  The  excitation  current  is 
uniform  across  the  length  and  width  of  the  rectangle.  The  eddy  current  describes  how  the  current 
is  disturbed  by  the  presence  of  a  magnetic  field. 

J tot  =  J ex  +  J eddy  (^.1) 

Applying  equation  (3.1)  to  the  case  of  uniform  excitation  current  and  an  eddy  current  density  that 
is  linear  in  x  results  in  the  illustration  of  Figure  9.  The  (+)  sign  indicates  current  is  flowing  into 
the  paper  and  a  (-)  sign  indicates  current  is  flowing  out  of  the  paper. 
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Figure  9  -  Cross-section  of  a  trace  illustrating  the  total  current  distribution 
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The  coordinate  system  for  a  rectangle  is  illustrated  by  Figure  10.  The  origin  is  located  at 
the  center  of  the  rectangle. 


Figure  10  -  Coordinate  system  for  the  top  view  of  a  rectangle 


The  author  in  [8]  gives  the  relationship  between  eddy  current  density  and  the  z- 
component  of  the  magnetic  flux,  using  Maxwell’s  equations,  as 

J eddy  =  °By  =  B:  (*  ')dx  ’  (3  -2) 

Note  that  the  x-component  of  the  eddy  current  is  neglected  since  the  length  of  the  rectangle  is 
assumed  to  be  much  larger  than  the  width.  A  reasonable  approximation  in  (3.2)  is  to  assume  that 
Bz  is  constant,  then  Jeddy  is  linear  in  x.  A  more  accurate  fonnulation  would  be  to  allow  J eddy  to 


vary  as  xa ,  where  a  is  a  constant,  determined  empirically,  that  is  dependent  on  the  technology 
being  used.  This  allows  equation  (3.2)  to  be  re-written  in  a  more  general  form  as 

deddy  =  jOJ<jB_Xay.  (3.3) 


The  author  in  [8]  uses  an  empirical  approximation  for  B~.  However,  the  magnetic  flux 
acting  on  a  given  rectangle  can  be  approximated  using  the  partial  inductance  matrix  found  from 
equation  (2.18).  The  total  flux  acting  on  rectangle  i  is 

N 

Vl='LMlJIj-  (3-4) 

j= i 

By  assuming  uniform  flux  distribution  over  the  area  of  rectangle  i,  the  magnitude  of  the  magnetic 
flux  density  can  be  approximated  as 


i  v 


B, 


7=1 


Wl 


(3.5) 
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The  components  of  the  partial  inductance  matrix  that  will  contribute  to  the  z  component  of  B, 
Mzi,  can  be  determined  by  the  right  hand  rule.  Orthogonal  rectangles  have  zero  mutual 
inductance.  We  neglect  the  self-inductance  term,  because  the  magnetic  field  generated  by 
rectangle  i  will  have  a  zero  average  z  component  inside  its  volume.  From  [11],  a  general 
expression  for  Mzi  is  given  as 

N 

M=,=  Z(-1  f«..  (36) 

n=l,n*i 


where  ain  is  determined  by  the  orientation  of  rectangle  i  relative  to  rectangle  n  and  is  equal  to 
zero  or  one.  Our  expression  for  Bz  on  rectangle  i  is 


B  -  MJ- 


WJ, 


(3.7) 


3.2  The  Skin  Effect 


In  general,  a  conductor  with  a  rectangular  cross  section  will  have  current  crowding  at  all 
of  the  surfaces  as  illustrated  by  Figure  11.  The  current  crowding  at  high  frequencies  is  due  to  the 
internal  inductance  being  higher  in  the  center  of  the  conductor  as  a  result  of  the  internal  magnetic 
field  being  largest  at  the  center.  The  current  tends  to  flow  through  the  path  of  least  impedance, 
which  is  the  path  of  least  inductance  at  high  frequencies. 


Figure  11  -  Cross-section  of  current  distribution 


An  exact  calculation,  in  terms  of  number  of  skin  depths  A  ( 1<  A  <2)  of  the  thickness  that 
the  current  in  the  conductor  penetrates  is  not  a  trivial  task  and  its  result  will  vary  from 
technology  to  technology,  depending  greatly  on  the  dielectric  properties  of  the  multilayer 
substrate  being  used.  Using  A  as  a  variable  for  the  number  of  skin  depths  the  current  takes  up,  the 
volume  current  density  in  the  rectangle  can  be  expressed  as 

Zj—Z 

J(x,y,z)  =  A*j(x,y)esy.  (3.8) 
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Here,  z„  is  the  distance  from  the  top  of  the  substrate  to  rectangle  i.  z  is  equal  to  zero  at  the  top  of 
the  substrate  and  is  negative  inside  the  substrate.  Integrating  equation  (3.8)  over  the  thickness  of 
the  conductor  gives 


'  j(x,y,z)dz  =  A* yj (x,  v) 


( 

S-Se5 

V  ) 


(3.9) 


The  effective  thickness,  /).//,  is  then  defined  as 


teff  =A*S 


f  -A 

l-e 5 


(3.10) 


V  J 

As  we  will  see,  teff  plays  a  critical  role  in  determining  the  frequency  dependent  resistance  and 
therefore  the  peak  Q  of  the  device.  This  allows  A  to  be  treated  as  a  degree  of  freedom  that  will  be 
determined  empirically  for  each  technology.  Using  Q  obtained  by  another  EM  solver  or  from 
measured  data  as  a  benchmark,  A  is  adjusted  so  that  the  peak  Q  obtained  by  this  model  matches 
that  of  the  EM  solver  or  measured  data.  For  the  technology  used  in  this  model,  a  value  of  A  =  1.1 
is  used. 


3.3  Calculating  the  Frequency  Dependent  Resistance 


The  total  power  dissipated  in  the  rectangle  can  be  expressed  in  terms  of  the  total  volume 
current  density  in  the  rectangle. 

(3-11) 

Jv  2 cr 1  1 

We  first  substitute  equation  (3.7)  into  equation  (3.3)  for  Bz,  plugging  that  result  into  equation 
(3.1)  and  equation  (3.11).  Equation  (3.11)  is  then  integrated  which  gives  the  total  power 
dissipated  as 


PJis,J^L\jj\W:  +  -  K 


2(7 


2a  +  1 


f  w  \2a+l  ( 


W, 

V  2  J 


Wt 

V  2; 


(3.12) 


The  variable  K  is  defined  as 


K  = 


(3.13) 


The  effective  resistance  of  the  rectangle  relates  the  amount  of  power  dissipated  in  the  rectangle 
to  the  amount  of  current  applied  to  the  rectangle. 
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2  P 

n  _  1  diss 

Keff 


2P 


diss 


14 


r 

[  (j  *dS  ) 

J  r,  y  ex  cross  J 

‘-'cross 

(3.14) 


Evaluating  the  integral  gives 


t  z  W 

L  1a  4*  *  ydxdz = 


(3.15) 


The  result  of  equation  (3.15)  and  equation  (3.12)  is  then  combined  into  equation  (3.14),  which 
results  in  Reff,md  for  a  given  rectangle  i  of  an  inductor  being 

^2a+l 

2"; 


R 


l 


eff  ,ind 


aWXf 


W,+- 


K 


2a  +  1 


fW-.  V“+1  f  W, 


v  ^  j 


(3.16) 


As  with  teff,  a  is  also  determined  empirically  by  the  comparison  of  Q  between  this  model 
and  another  reliable  benchmark.  The  procedure  for  determining  a  is  to  start  off  with  a  =  1 .  The 
value  for  te/f  is  varied  first,  as  its  value  will  affect  Q  the  most.  The  value  for  teff  will  always  be 
between  one  and  two.  If  adjusting  %-does  not  give  a  correct  value  for  peak  Q,  tejjw\\\  be  set  to 
the  value  that  gives  the  best  Q.  Next,  the  value  for  a ,  is  adjusted.  Its  value  should  not  be  too  far 
from  1,  as  the  linear  approximation  for  the  eddy  current  distribution  is  fairly  accurate.  For  this 
model,  the  linear  approximation  for  the  eddy  current  distribution,  a  =  1,  is  used. 

The  rectangles  of  the  inductor  will  experience  an  internal  inductance  because  of  the 
magnetic  flux  internal  to  the  conductor.  Paul  [3]  gives  the  per-unit-length  internal  inductance  of 
a  circular  wire  as 


L  =  — 


2^/0 


v8  n  j 


(3.17) 


for  high  frequencies  (rw»  S)  and 


rillt  _  Mo 

8ft 


(3.18) 


for  low  frequencies  (rw«S).  Here,  rw  is  the  radius  of  the  wire.  28  is  recognized  as  the  effective 
thickness  of  the  wire.  As  the  conductor  is  rectangular,  it  is  unknown  what  should  be  substituted 
into  rw.  However,  if  we  assume  that  at  DC,  the  per-unit-length  internal  inductance  should  be 
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equation  (3.18),  then  it  is  reasonable  to  approximate  the  per-unit-length  internal  inductance  for 
all  frequencies  as 


jint  _  leff 


t,  ' 

t  K 


Mo. 

8;r 


(3.19) 


because  at  low  frequencies,  tejgr  is  approximately  equal  to  t.  The  total  reactance  of  the  rectangle 
can  be  calculated  as 

Xeff=a>Lmtl.  (3.20) 

Chapter  4  Deriving  the  Inductor  Lumped  Pi-Model 

The  techniques  of  Chapters  2  can  be  used  to  construct  a  PEEC  distributed  model  for  any 
planar  monolithic  inductor.  Shown  in  Figure  12  is  a  flowchart  outlining  the  progression  from 
Maxwell’s  equations  to  an  inductor  lumped  pi-model. 


Figure  12  -  Inductor  lumped  pi-model  derivation  flowchart 
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In  this  chapter,  we  describe  a  method  to  simplify  the  distributed  model  down  to  a  lumped 
pi-model  for  a  spiral  inductor.  A  lumped  pi-model  makes  it  easy  for  the  designer  to  visualize  the 
parasitics  inherent  to  the  spiral  inductor.  The  inductor  lumped  pi-model  we  use  is  commonly 
used  for  integrated  inductors  and  is  shown  in  Figure  13.  The  loss  in  the  substrate  due  to  the 
electrostatic  coupling  between  metal  and  ground  is  captured  by  Rsi  and  RS2.  The  capacitance 
caused  by  this  electrostatic  coupling  through  the  substrate  is  captured  by  Csi  and  CS2.  The 
capacitive  coupling  due  to  the  lossless  oxide  is  represented  by  Coxi  and  C0X2-  The  lumped  model 
used  by  [11]-[13]  does  not  account  for  this  oxide  capacitance  separately  from  Csi,  Rsi,  CS2,  and 
RS2.  By  taking  into  account  the  effect  due  to  Coxi  and  C0X2  separately  from  the  rest  of  the 
substrate  admittance,  we  are  able  to  accurately  model  the  device  at  lower  frequencies.  If  Coxi  and 
COX2  are  absent,  we  effectively  have  a  low  impedance  to  ground  through  Rsi  and  RS2  at  low 
frequencies,  which  if  physically  not  accurate  since  the  lossless  oxide  would  make  the  impedance 
very  large. 

As  mentioned  in  Chapter  1 ,  there  are  approximate  empirical  and  semi-analytical  methods 
([1],  [6],  [8],  [15],  [16])  used  to  calculate  the  component  values  of  the  lumped  pi-model.  As 
illustrated  in  the  flowchart  in  Figure  12,  we  use  two  low  frequency  approximations  to  the 
distributed  model  similar  to  that  of  [11]-[13]  to  calculate  Rs,  L,  Csi,2  and  Rsi,2-  We  use  a  high 
frequency  approximation  that  uses  the  overlap  capacitance  of  the  inductor  bridge  to  approximate 
a  value  for  C0  and  Rq. 


Figure  13  -  Inductor  lumped  pi-model 


Approved  for  public  release;  distribution  is  unlimited. 

18 


ASCSII-2  contract  #F33615-00-D-1726-0002 


CDRL  #A001  Final  Report 


4.1  First  Low  Frequency  Approximation 

The  first  low  frequency  approximation  involves  recognizing  that  at  low  frequencies,  the 
complex  capacitors  are  approximately  open  circuits.  By  applying  this  approximation  to  the 
PEEC  circuit  of  Figure  8,  the  circuit  of  Figure  14  results. 


Figure  14  -  Distributed  circuit  following  a  low  frequency  approximation 


For  a  circuit  with  N  series  rectangles,  it  is  easily  observed  that  the  total  series  impedance  of  this 
circuit  is 


The  form  of  equation  (4.1)  is 


N  N  N 


i= 1  i= 1  y=l 


-  R  ,  +  sL 


which  is  the  same  form  as  the  branch  of  the  pi-model  shown  in  Figure  15. 


(4.1) 


(4.2) 


Figure  15  -  Pi-circuit  neglecting  capacitances 


Thus  the  components  of  Figure  14  can  be  found  as 

N 

R  =YR . 

s  /  J  i 

i= 1 
N  N 

i= 1  7=1 


(4.3) 

(4.4) 
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4.2  Second  Low  Frequency  Approximation 

The  second  low  frequency  approximation  does  not  neglect  all  of  the  capacitors.  Because 
the  impedance  of  an  inductor  is  much  smaller  than  the  impedance  of  a  capacitor  at  low 
frequencies,  the  complex  capacitances  in  parallel  with  a  series  RL  combination  are  neglected.  All 
of  the  complex  capacitances  connected  to  ground  remain.  When  this  approximation  is  applied  to 
the  PEEC  circuit  of  Figure  8,  the  result  is  the  circuit  of  Figure  16.  The  currents  through  each  of 
the  rectangles  of  the  inductor  are  assumed  to  be  approximately  equal  because  of  the  negligible 
current  in  the  substrate.  This  allows  us  to  lump  the  mutual  inductances  with  the  self-inductances. 


Port  2 


Figure  16  -  PEEC  circuit  after  second  low  frequency  approximation 


4  =!>,;/  (4-5) 

7=1 

The  next  step  is  to  short  port  2  and  calculate  the  input  admittance  Yu.  We  ignore  terms 
with  s  to  the  power  of  three  or  greater  since  we  are  using  a  low  frequency  approximation.  This 
gives 


E, 


Cm,CA 

kCox1+C's1 


■  +  CC^ 


c  c 

'~'ox2'~'s2 

C  +C' 

^ox2  T 


+  a , 


cm3c:3 

Cox3+C's3J 


+  - 


(4.6) 


where  Zseries  is  found  from  equation  (4.1)  and 


N 


series 


slLiL 

7=1 

L 


(4.7) 
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if  we  assume  sLt  »  Rj  as  in  [11]  to  make  the  substrate  impedance  independent  of  frequency. 
The  form  of  equation  (4.6)  matches  the  form  for  the  input  admittance  looking  into  the  circuit  of 
Figure  17. 


Figure  17  -  Low  frequency  pi-circuit  with  port  2  shorted 


Y  =  sY'  + 

-'ll  J  wi  T 


R  +  sL 


Equating  equation  (4.6)  to  equation  (4.8)  for  N  rectangles  gives  the  substrate  admittance  Y'sX . 
Since  s=jco  and  keeping  in  mind  that  C'si  is  complex, 


N  N 

y;t=2aX  =Ea.7rjT7ir=i®  M?/,)  +>h*Im(f.',)} 

*=1  *=1  ^oxi^^si  (4.9) 

=  >{Re(rl)}-«{im(r;l)} 

where  Y'sX  is  the  total  complex  admittance  to  ground  at  port  1  as  illustrated  in  Figure  17.  Y'si  is  the 
total  complex  capacitance  to  ground  of  rectangle  i  and  is  found  from  the  capacitance  matrix  as 

N 

K=Xcv’  (4-10) 

7=1 

where  c,y  represents  the  individual  terms  of  [C]. 

In  a  manner  similar  to  [15]  and  [16],  we  use  a  parallel  plate  approximation  for  the  oxide 
capacitance.  This  gives  the  oxide  capacitance  in  terms  of  the  area,  thickness,  and  permittivity  of 
rectangle  i  as 

coxi  =  AiyL>  (441) 

^OX 

where  At  is  the  surface  area  (neglecting  thickness)  of  rectangle  i.  We  make  the  assumption  that 
the  oxide  capacitance  is  equal  at  both  ports;  therefore,  we  can  take  the  sum  of  the  oxide 
capacitance  due  to  all  the  rectangles  and  divide  it  by  two  to  calculate  the  oxide  capacitance  at 
port  1  and  port  2. 
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C  ,  =C  , 

ox, pi  ox,p2 


i  N 

=  -fc  . 

^  /  J  OX! 


2tr 


(4.12) 


We  can  express  F'  as  the  series  combination  of  the  admittances  due  to  C'sM  and  C  , 
which  are  defined  in  Figure  17. 

Q,A,,(Qj+c„,)+c„X.|  I  c£*c„ 


-<U< 


(c^  +  c^+c^ 


(4.13) 


where  C71  =  Im  {CsM }  and  CR ,  =  Rc  { C’s'(y/j ,  j  .  Equating  the  real  and  imaginary  parts  of  equation 

(4.9)  with  the  real  and  imaginary  parts  of  equation  (4.13)  gives  us  two  equations  to  solve  for  the 
two  unknowns,  Qi  and  CR  \.  From  Qi  and  CRi \  we  use  the  complex  capacitance  technique  to 
extract  the  component  values,  Csub\  and  Rsub\- 

(4.14) 

1 


c  =c 

'-'sub  1 


R sub\ 


coC 


(4.15) 


i,  l 


This  technique  will  work  for  any  co  >  0.  A  similar  procedure  is  applied  to  find  the  component 
values  at  port  2. 

N  N 

%  =  £  A  n  =£  A  xAtfr  =  j0>  { Re  ( f;2 )}  +  jo,  {j  » Im  ( )} 

i= 1  i=l  Co.«  +  (4.16) 

=  j(0  j  Rc  ( T  o  )|  - (!)  j  Im  ( F  x ) 


£  = 


y=i 


L 


c  =  c 

^  sub!  '-'R,  2 


R subl 


coC 


(4.17) 

(4.18) 

(4.19) 


1,2 


4.3  High  Frequency  Approximation 

For  an  accurate  calculation  of  the  self-resonant  frequency,  the  most  important  component 
of  the  equivalent  pi-model  apart  from  the  inductance  is  the  interwinding  capacitance  C'0 ,  the 

parallel  combination  of  R0  and  C0  in  Figure  13.  While  the  frequency  dependent  resistance  plays 
an  important  role  in  determining  the  value  of  the  peak  Q  of  the  inductor,  the  interwinding 
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capacitance  determines  when  the  inductor  Q  becomes  zero  and  the  inductor  behaves  as  a 
capacitor  rather  than  an  inductor.  We  replace  the  complete  distributed  model  for  three 
consecutive  rectangles  with  the  simplified  version  shown  in  Figure  18  for  4  rectangles,  where  the 
inductors  are  treated  as  open  circuits.  Note  we  also  replace  the  admittances  to  ground  from  each 
rectangle  with  the  lumped  values  at  the  ports  found  with  the  second  low  frequency 
approximation. 


Figure  18  -  High  frequency  approximation  for  four  series  rectangles 

Our  method  is  based  on  a  common  technique  used  in  other  analytical  pi-models,  such  as 
[15]  and  [16],  which  use  the  overlap  capacitance  of  the  bridge,  illustrated  in  Figure  19,  to 
determine  the  interwinding  capacitance. 


Figure  19  -  Illustration  of  bridge  overlap  capacitance  -  2.5  turn  inductor 
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The  methods  in  [15]  and  [16]  use  a  parallel  plate  approximation,  whereas  we  use  the  [C] 
matrix,  which  is  more  accurate  in  that  it  accounts  for  fringing.  In  order  to  justify  the  method  of 
using  the  overlap  capacitance  to  detennine  the  total  interwinding  capacitance,  we  must  show  that 
the  portion  of  the  inductor  that  overlaps  the  bridge  is  close  to  the  port  1  voltage  and  that  the 
bridge  is  close  to  the  port  2  voltage  at  resonance.  Once  this  is  accomplished,  we  can  ignore  all  of 
the  turn-to-turn  interwinding  capacitances  because  most  of  the  turns  are  at  the  same  potential  and 
the  overlap  capacitance  will  dominate.  After  inspecting  the  voltage  along  the  inductor  for  various 
geometries  by  solving  the  above  equations,  we  consistently  find  that  the  voltage  along  the  first 
half  of  the  inductor  rectangles  remained  fairly  constant  and  that  most  of  the  voltage  drop 
occurred  along  the  inner-most  turns  of  the  inductor  after  the  last  overlap  of  the  bridge.  The  plot 
of  the  voltage  along  the  inductor  at  resonance  is  shown  below  for  a  1.5  turn  and  3.5  turn  inductor 
in  Figure  20  and  Figure  21,  respectively,  and  it  shows  that  the  overlap  rectangles  are  close  to  the 
same  voltage  as  port  1.  Similar  results  were  found  for  other  N-turn  inductors. 


1  M 

i  i 


■■-rUu-L  l-tlt  v«  IlhJilCliV 


2 1 


Figure  20  -  Inductor  voltage  vs.  Inductor  rectangle  for  1.5  turn  inductor 

hMiuniU  VfiiFl  U|>K  arm  lllilaii'  Enr 


Figure  21  -  Inductor  voltage  vs.  Inductor  rectangle  for  3.5  turn  inductor 
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As  a  result  of  these  findings,  we  can  express  the  total  interwinding  capacitance  as  the 
parallel  combination  of  each  bridge  overlap  capacitance.  A  general  expression  for  the 
discretization  scheme  shown  in  Figure  3  and  used  in  this  model  is  given  as 

7-9 

co  =  X  (c(fi  r) + c{i + 1,  r)}  (4.20) 

i=6 

stepS 

y  =  8A  +  3  (4.21) 

c(i,j)  =  -C(iJ)  (4.22) 

with  N  being  the  number  of  turns  of  the  inductor. 

We  use  the  complex  capacitance  technique  to  extract  the  real  capacitance  and  resistance 
values  from  the  complex  capacitance. 

C0=reaI(C'0)  (4.23) 


oj(imag(C'l )) 


(4.24) 


Chapter  5  Comparison  with  Momentum  and  Measured  Data 
5,1  Inductor  Results 

A  set  of  five  spiral  inductors  was  designed,  with  the  1 .5  turn  design  being  done  by  Idstein 
[11],  with  turns  varying  from  1.5  to  5.5  turns  as  described  by  Table  1.  The  substrate  used  is 
based  on  a  0.18  um  FDSOI  CMOS  technology,  a  process  utilizing  a  high  resistivity  (2000  Q- 
cm.)  bulk  silicon  substrate.  A  description  of  the  layer  stack  and  its  properties  is  given  in  Tables  2 
and  3. 


Table  1  -  Inductor  geometry  summary 


Inductor 

Inductance 

Number  of 

Outer  diameter 

Trace  width 

Trace  spacing 

# 

(nH) 

turns 

(um) 

(um) 

(um) 

1 

1 

1.5 

248 

28 

4 

2 

2.5 

2.5 

248 

14 

4 

3 

1.6 

3.5 

146 

10 

4 

4 

3 

4.5 

190 

10 

5 

5 

4.9 

5.5 

210 

10 

4 
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Material 

Thickness  (um) 

Dielectric  Constant 

Conductivity  (S/m) 

Oxide 

7.75 

4 

0 

Silicon 

675 

11.7 

0.5 

Table  2  -  Substrate  dielectric  properties 


Metal  Layer 

Thickness  (um) 

Distance  from 

metal  center  to 

top(um) 

Conductivity 

(S/m) 

Metal  1 

.6 

6.15 

2.08E+07 

Metal  2 

.6 

4.55 

2.08E+07 

Metal  3 

2 

1.75 

2.64E+07 

Table  3  -  Substrate  metal  properties 


Using  the  concepts  discussed  in  this  report,  the  component  values  of  the  lumped  pi-circuit 
of  Figure  13  were  calculated.  The  lumped  pi-circuit  was  then  simulated  using  the  Agilent  ADS 
schematic  tool.  The  ADS  schematic  showing  the  pi-circuit  and  its  component  values  for  inductor 
#1  is  given  in  Figure  22.  The  box  in  the  middle  of  the  schematic  titled  “SNP1”  contain  the  2-port 
S-parameters  of  the  frequency  dependent  impedance.  The  2-port  S-Parameters  were  computed 
from  100  MHz  to  the  first  resonant  frequency  of  each  inductor.  The  inductors  were  also 
simulated  using  ADS  RF  Momentum,  which  is  a  Method  of  Moments  (MoM)  EM  simulator. 
Measured  data  was  obtained  for  inductor  #1,  the  1.5  turn  spiral.  The  layout  of  the  test  structure  is 
illustrated  by  Figure  23.  The  pads  are  arranged  in  a  ground-signal-ground  configuration  and  2- 
port  S-Parameter  measurements  from  1-20  GHz  were  performed  using  on-chip  cascade  probes. 
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Figure  22  -  ADS  schematic  for  inductor  #1 


Figure  23  -  1.5  turn  inductor  (a)  and  open  (b)  test  structure 


The  measured  data  was  de-embedded  by  subtracting  the  Y-Parameters  (Yu  +  Y12,  Y22  + 
Y12)  of  the  de-embed  open  (Figure  23b)  from  the  measured  Y-Parameters.  This  effectively 
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removes  the  parasitic  pad  capacitance  from  the  measurements.  For  the  other  inductors,  the  2-port 
S-Parameters  generated  from  the  Momentum  simulation  were  used  as  the  benchmark  to  compare 
our  model  to. 

The  differential  impedance,  which  describes  the  series  impedance  connecting  ports  1  and 
2,  is  defined  using  the  Z-Parameters  of  the  2-port  inductor. 

Zdiff  =  Z\\  +  z22  —  Z12  —  Z21  (5.1) 

Using  the  differential  impedance,  the  effective  inductance  of  the  inductor  can  be  computed  at 
frequency  /  as 

imag(Zdiff) 


Leff  -  ' 


2  ttf 


(5.2) 


The  inductor  Q  can  also  be  computed  using  the  differential  impedance  as 

imag(ZdW) 


Q  = 


real(zdiff) 


(5.3) 


The  values  for  Q  and  Leff,  as  well  as  the  2-port  S-Parameters  generated  by  Momentum 
and  our  model  were  then  compared.  Again,  we  were  able  to  compare  to  measured  data  for 
inductor  #1.  As  can  be  seen  from  Figure  24  -  Figure  28,  reasonable  agreement  between 
Momentum  and  our  model  is  achieved.  Additionally,  our  model  agrees  well  with  the  de- 
embedded  measured  data  in  predicting  the  frequency  and  value  of  the  peak  Q ,  as  well  as 
accurately  modeling  effective  inductance  and  S-parameters,  as  can  be  seen  in  Figure  29. 
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.lo- 


Figure  24  -  Simulated  results  inductor  #1 


Si  1  -  25  lurn  inductor  S12  -  2.5  turn  inductor 


Figure  25  -  Simulated  results  inductor  #2 
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Figure  26  -  Simulated  results  inductor  #3 


St  I  •  4.5  turn  inductor  S12  •  4.5  turn  inductor 


Figure  27  -  Simulated  results  inductor  #4 
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Si  1  •  5.5  turn  inductor  S12  •  5.5  turn  inductor 


Figure  28  -  Simulated  results  inductor  #5 


811  S12 


Figure  29  -  De-embedded  measured  results  for  inductor  #1 
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The  biggest  discrepancy  observed  is  in  the  prediction  of  Q  at  higher  frequencies.  It  is 
typical  for  Q  to  vary  slightly  from  one  model  to  the  next,  as  the  methods  of  capturing  the  ohmic 
loss  in  the  metal  are  usually  different.  The  most  important  aspect  of  Q  is  the  magnitude  and 
frequency  of  its  maximum,  since  this  is  typically  in  the  frequency  range  that  the  device  will  be 
used.  In  comparing  our  model  with  Momentum,  the  peak  Q  is  reasonably  close,  and  the  peak  Q 
frequencies  are  nearly  identical  which  leads  to  a  high  degree  of  confidence  in  both  models.  Slight 
deviations  in  the  effective  inductance  values  are  most  likely  due  to  approximations  made  in  the 
modeling  of  internal  inductance.  The  slight  over-estimation  in  the  resonant  frequency  is  most 
likely  due  to  the  method  in  detennining  the  interwinding  capacitance,  which  neglects  some  of  the 
turn-to-turn  capacitances  and  relies  only  on  the  bridge  overlap  capacitance  for  its  value. 

Chapter  6  Summary  and  Conclusions 

In  this  report,  the  technique  proposed  in  [11]-[13]  for  modeling  spiral  inductors  with  a 
simple  lumped  pi-circuit  was  improved  upon.  In  addition,  the  bandwidth  of  this  model  goes  up  to 
the  first  resonant  frequency,  a  significant  improvement  in  the  bandwidth  achieved  previously. 

As  in  [11]-[13],  the  PEEC  technique  was  used  to  construct  the  distributed  circuit,  upon 
which  the  lumped  pi-circuits  were  based.  A  low  frequency  approximation  similar  to  that  in  [11]- 
[13]  was  applied  to  the  distributed  circuit  to  arrive  at  values  for  the  substrate  impedance  and 
series  impedance  of  the  inductor.  The  method  in  this  work  differs  slightly  in  that  it  accounts  for 
the  lossless  oxide  capacitance.  At  higher  frequencies,  a  different  method  was  used  to  determine 
the  interwinding  capacitance.  This  method  has  also  been  used  in  [15]  and  [16],  but  we  justify  its 
application  by  inspecting  the  voltage  along  the  inductor  near  resonance.  As  previously 
mentioned,  the  interwinding  capacitance  must  be  calculated  accurately  in  order  to  correctly 
predict  the  resonant  frequency. 

In  calculating  the  partial  capacitance  matrix,  [11]-[13]  used  numerical  integration  which 
proved  to  be  very  time-consuming.  A  method  employing  FFTs  as  discussed  by  [5]  was  used  in 
this  work  to  speed  up  the  calculation  of  [C]. 

As  demonstrated  in  Chapter  5,  good  agreement  in  Q,  Leff,  and  2-port  S-parameters  up  to 
the  first  resonance  is  obtained  with  the  inductor  model  presented  in  this  work  and  ADS  RF 
Momentum.  We  also  verify  the  accuracy  of  the  model  by  comparing  to  the  de-embedded 
measurements  for  a  1.5  turn  1  nH  inductor. 
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A  big  improvement  that  could  be  made  to  this  model  would  be  to  account  for  eddy 
current  losses  in  the  substrate  when  calculating  the  partial  inductance  matrix  \M\.  To  accomplish 
this,  a  multi-layer  Gm  would  need  to  be  calculated,  as  mentioned  in  Chapter  2.  The  technique 
outlined  in  [19],  using  image  theory,  could  be  used  as  a  guideline  for  calculating  Gm.  One 
potential  problem  with  using  a  Gm  that  accounts  for  loss  is  that  the  matrix  \M]  will  become 
frequency  dependent,  which  leads  to  a  somewhat  more  complicated  model. 

Another  improvement  to  the  inductor  model  would  be  to  allow  for  more  general 
geometries,  such  as  octagonal  spirals  or  circular  spirals.  As  the  shape  of  a  spiral  inductor 
approaches  circular,  the  efficiency  in  terms  of  inductance  per-unit-area  increases,  leading  to 
more  efficient  use  of  real-estate  in  the  integrated  circuit  design.  The  only  change  that  needs  to  be 
done  is  to  modify  the  existing  computer  program  (MATLAB). 

The  model  discussed  in  this  report  is  based  on  rigorously  derived  equations  that  are  true 
to  the  physics  of  the  device.  The  S-parameter  data  generated  in  the  ADS  schematic  can  be 
incorporated  into  other  circuit  simulators  such  as  Cadence  and  used  in  the  design  of  complex 
integrated  circuits  such  as  LNAs,  matching  networks,  LC  tanks  for  oscillators,  etc.  Since  this 
model  can  be  run  in  a  very  short  period  of  time,  the  design  cycle  time  can  be  greatly  reduced 
compared  to  the  time  involved  with  other  numerical  EM  solvers. 
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Appendix  A:  A  General  Approach  to  Solving  the  Electrostatic  Green’s  Function  for  a 
Multilayered  Substrate 


Z  Free  Space 


PEC  Ground 


source  location 


Figure  30  -  Profile  of  F-layered  substrate 


With  reference  to  [5],  [10],  and  [18],  we  consider  the  derivation  of  the  multilayer 
electrostatic  Green’s  function,  Ge,  for  the  L-layered  structure  in  Figure  30.  The  dimensions  of  the 
structure  shown  in  Figure  30  are  a  in  the  x-direction  and  b  in  the  y-direction.  Ge  (/•;/'’)  is  found 

by  inverting  Poisson’s  equation  and  gives  the  potential  at  the  observation  point  r  due  to  a  point 
charge  source  placed  at  a  point  r' .  In  our  case,  all  metal  layers  are  in  the  top  layer,  which  means 
both  the  source  and  observation  points  will  be  in  the  top  layer.  However,  for  generality  we  will 
perform  the  derivation  assuming  the  source  location  can  be  in  any  layer. 

Poisson’s  equation  in  equation  (A.l)  can  be  written  in  the  Green’s  function  fonn  of 
equation  (A. 2). 

VV  =  -—  (A.l) 

£ 

/  r\  S(x-x')S(v-v')S(z-z') 

V2Ge(r;r')  = - - - — — (A.2) 
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where  the  subscript  i  indicates  the  source  layer.  The  form  of  Ge  is 

Ge  =X(x;x')Y(y;y')Z'(z;z').  (A3) 

We  substitute  A  Y,  and  Z'  into  equation  (A.2)  which  gives 

d2  d2  d2  S(x  —  x')S(y—y')S(z  —  z') 

YZ — -X  +  ZX — -Y  +  XY — -Z  = - - - ’  Ky  y’  y (A.4) 

dx  dy  dz~  sj 

The  boundary  condition  on  the  side-walls  of  the  structure  in  Figure  30  is  h»E  =  0 .  To  satisfy 


J717T  TUT 

this,  we  let  X  =  cos  - x  and  Y  =  cos  —  v  .  Plugging  in  the  functions  for  X  and  Y  we  re- 

v  a  J  l  b  ' 


write  equation  (A.4)  as 


vn  filJL  HJL  U  . 

2^2^iC0S  — x  cos  — y  x i — ~ 

m=0  n=0  V  a  )  \b  )  dz 

S{x-  x')8  (y  -  y')8{z  -z') 


7-  (m^X  frtiTt 


Z'  = 


YY17T  TUT 

We  multiply  both  sides  by  cos  - x  cos  —  v  and  integrate  both  sides  over  x  and  v  from  0 

{  a  J  {  b  J 


to  a  and  0  to  b,  respectively,  and  this  gives 


ab  d2  r7l  f  V  f  nxn  \  \r7,  \ 

— ) — tZ  -  1 —  ]  I Z  }  =  - 

4  dz  \  a  )  \  b 


S(z-z')  ( nm  ( rm  ? 

- COS  - X  COS  - V  (A.6) 

s,  [a  J  l  b  '  J 


TY17T  TUT 

We  substitute  Z'  =  Zmn(z;z') cos  - x'  cos  — y  into  equation  (A.6)  to  simplify  to 

\  a  )  \  b  "  ) 


^  d^_f2  7  S(z-z') 

.  j  2  hnn  ^ mn  \Z’ Z  ) 

4  dz  £: 


,2  f  mxX  ( nn  V 


where  t2mn  =  -  +  —  .  From  equation  (A. 7)  we  can  solve  for  Zmn  (z;zr)  using  the  method 

{  a  J  {  b  J 


of  solving  a  one-dimensional  Green’s  functions  outlined  in  [18].  The  method  is  summarized 


below. 


A  solution  to  equation  (A. 7)  has  the  form 


zmn(Gz')  =  Cke^  +  Dke~,-Z 


in  layer  k.  A  recursive  relation  is  obtained  between  the  coefficients  in  adjacent  layers  using  the 


boundary  conditions.  At  the  boundary  between  each  interface,  the  voltage  and  nxD  is 
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continuous.  In  addition,  the  voltage  at  the  bottom  is  zero  due  to  the  ground  plane.  At  the  top,  we 
have  zero  normals  field  due  to  the  magnetic  wall  condition.  This  recursive  relation  can  be 
represented  by  a  matrix  [A], 


C, 


*+i 


D 


k+ 1. 


=  [A\ 


Q 

D, 


(A-9) 


where  [A]  can  be  found  from  the  boundary  conditions  as 


M= 


f  e  ' 

1  +  -^ 

V  ek+ 1 J 


1  — — A 


V  £k+ 1 J 


-2  tm,A 


1 

(  e  > 

1 

2  tmndk  i 

e  k  — 

1+  £k 

2 

V  £k+ 1  ) 

2 

V  G+l  y 

(A.  10) 


Next,  we  introduce  the  self-adjoint  operator  L. 

i=TL(z)A 

dz  dx 


-q(x) 


We  want  a  Green’s  function  that  satisfies  an  equation  of  the  form 

[L  +  Aw(z)]Zmn  (z;  z')  =  -8(z  -  z') 
together  with  the  boundary  conditions 

a  ,  p(z)Z'  (z:z')+  B  ,Z  (z;z')\  ,  =  0, 

'  mn  V  ?  /  A  a.b  mn  V  ’  '\z=a.b  ’ 


(A.  11) 


(A.  12) 


(A- 13) 


where  L  is  the  self-adjoint  operator  of  equation  (A.l  1).  The  solution  to  equation  (A.  12)  is  divided 
into  two  intervals. 

U(z)T(z') 


Z  (z:z')  = 

mn  \  ’  / 


C 

U(z')T(z) 

C 


z<z 

z>z 


(A.  14) 


where  C  is  a  constant  and  U  and  T  are  two  independent  solutions  of 

U  (z) 

[L  +  Aw(z)\  )  =0 

T(z) 


(A.  15) 


with 


aaP(z)U'(z)  +  pcfj  (z)  =  0  [  z= 
ahP(z)T'(z)  +  fibT  (z)  =  0 1  z=b 


(A.  16) 


and 
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C  =  p{z')[T(z')U\z')  -  T'(z')U(z')\ .  (A.  17) 

Equation  (A.  16)  is  necessary  to  ensure  that  equation  (A.  13)  is  satisfied.  To  use  this  method  for 
equation  (A. 7),  we  re-write  it  as 


where 


dbS  d  2  ry  ,  f X 

- —  tt-L  Zm,Xz;z)  =  -S(z-z) 
A_„  dz 


4  m,n^0 

Am„  =  <  2  (m  =  0,  n  ^  0 )or(m  ^  0,  n  =  0) . 
1  m =n= 0 


(A.  18) 


(A.  19) 


To  put  equation  (A.  18)  into  the  same  form  as  equation  (A.  12),  we  recognize  that  p(z)  =  ■ 


w(z)  =  1 ,  q(z )  =  0 ,  and  A  =  -t2mn .  In  addition,  aa  =  0 ,  fda=  1 ,  =  0 ,  and  ah  =  — —  to  satisfy 


the  boundary  conditions  of  equation  (A.  16). 

We  proceed  with  solving  for  T(z)  and  U(z')  for  z>z.  Using  equation  (A.8)  in  matrix 


form,  we  solve  for  T(z)  as 


T(z)  =  \e ^  e * 


(A. 20) 


where  k  >  i .  We  can  find  Ck  and  D,  using  the  boundary  condition  at  z  =  d/_  in  conjunction  with 
the  recursive  relation,  applied  from  the  top  down  towards  the  k1h  layer.  From  the  boundary 


condition  at  dr  we  find 


Dl  =  CLel 


(A.21) 


We  can  set  DL  to  1  since  it  will  drop  out  in  a  ratio  with  D/  later  on.  This  allows  us  to  find 


as 

C"  =iA] -\AmT\..[AL-2X\Al^  2]  CL  (A.22) 

L Dk\  ie  -  J 

We  plug  the  result  for  equation  (A.22)  into  equation  (A.20)  to  find  T(z).  Next  we  need  to  find 
U ( z' )  for  z>z. 
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(A.23) 


The  superscript  /  in  the  coefficients  C  and  D  indicates  that  the  source  location  z  is  below 
the  observation  point  z  .  The  coefficients  C\  and  D.  for  U (z')  are  found  in  a  manner  similar  to 
how  the  coefficients  Ca  and  Da  were  found  for  T(z)  except  we  start  from  the  bottom  and  apply 
the  recursive  relation  towards  the  top.  From  the  boundary  condition  at  z  =  0  where  we  have  a 
ground  plane,  we  can  relate  the  coefficients  C;  and  D/  by  plugging  z  =  0  into  equation  (A. 8)  and 
setting  Zmn  (0;z')  =  0 .  From  this  we  get 


c1=-a 


(A. 24) 


We  can  set  D;  =  1  for  the  same  reason  we  can  set  Dl=  1  previously.  This  gives  us 


c; 

d! 


=  [Al_l][Ai_2\...[Ai] 


1 

-1 


C, 


(A. 25) 


The  values  for  the  coefficients  found  in  equation  (A.25)  are  then  plugged  into  equation  (A.23). 

Up  to  this  point,  half  of  the  solution  to  equation  (A.  14)  has  been  found.  The  next  step  is 
to  solve  for  T(z')  and  U(z )  for  the  case  z  <  z  .  For  this  case,  layer  k  <  layer  i.  U(z )  is  found  as 


r  tm„z  -tm„z  ~i 

— i 

_ i 

pk_ 

(A. 26) 

The  coefficients  Ca  and  Da  are  found  using  the  recursive  relation  applied  from  the  bottom  layer 
to  the  k‘h  layer,  with  the  aid  of  the  boundary  condition  Zmn  (0;z')  =  0 . 


re,- 

'  1 1 

Dk\ 

=  [4-i][A-2]-[A][4] 

-l 

ci 

(A.27) 

Again,  we  plug  the  result  of  equation  (A.27)  into  equation  (A. 26).  T(z')  is  found  for  z  <  z  as 


T(z')  =  [C 


(A.28) 


The  superscript  u  in  the  coefficients  C  and  D  indicates  that  the  source  location  z  is  above  the 
observation  point  z  .  The  coefficients  C“  and  D,"  are  found  by  applying  the  recursive  relation 


with  [A],  starting  again  from  the  z  =  cl/  ,  in  conjunction  with  the  boundary  condition  at  z  =  di  and 
working  downwards  towards  layer  i.  This  gives 
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=[Air[AMr...[AL-,r\AL_xr 


r 

2t„„d,  L  ' 


(A.29) 


Again,  we  plug  the  result  of  equation  (A.29)  into  equation  (A. 28). 

An  expression  for  C  in  equation  (A.  14)  is  found  using  equation  (A.  18).  However,  the 
expression  found  for  C  would  be  broke  into  two  intervals,  one  for  z  <  z'  and  one  for  z  >  z' .  We 
introduce  a  variable  C ,  as  well  as  coefficients  C, .  D  ,  C  ,  and  D  to  allow  us  to  write  Zmn(z\z  ) 


in  one  expression. 


Z  (z;z')  =  ir« 

mn  V  5  /  ^ 


:][c<]re'”z>  e^Ac> 

JL^<JL  \D>1 


=  Ge(z;z') 


(A. 30) 


where  z<  clearly  means  the  use  of  z  if  z  <  z  and  the  use  of  z  if  z  >  z' .  We  also  have 

c  \Ck  k<i,z<z' 

<  [C?  k>i,z>z' 

D  =\Dk  k-^z~z' 

<  I  D.  k>i,z>z' 

c  \Ck  k>i,z>z' 

>  [C”  k<i,z<z' 

D  ={Dk  k~kz~z' 

>  |D"  k  <i,z<  z' 


(A.31) 


(A.32) 


(A. 33) 


(A. 34) 


The  variable  C  is  defined  as 


0  1 


-Ul  -1  ][Z1f[Z2f...[4_2f[4_1f  ^  q  WTU+J1 


t^rWir1  2ti . 


(A. 35) 


We  tested  this  method  by  applying  it  to  the  substrate  under  consideration  in  this  report. 
For  this,  we  set  i=j= 2  and  L= 2  since  we  have  a  two-layer  substrate  with  both  metal  layers 
(source  and  observation  points)  in  the  top  layer.  Without  getting  into  the  details,  we  find 
Zm„(z;z')  as  follows.  For  z  <  z  , 


v  _  ( C'e'™z  +  D'e^'"z)(e'“z'  + 

A mn  (z’ z  )  ~  Onhr 

^  -D'2) 


(A.36) 
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Equation  (A.36)  matches  the  result  found  in  [10]  found  for  the  case  z<z  .  For  z  >  z' 

(C'e'"2'  +  D^™Z')(V"“Z 


Z  (z;z')  = 

mn  V  5  / 


labs ; 


(A. 37) 


A. 


'  t  (C'e 

mn  V  2 


/  ltd, 


-D‘) 


This  also  matches  the  result  found  in  [10]  for  the  case  z  >  z' . 

The  above  derivation  works  for  all  values  of  m  and  n  except  for  m  =  n  =  0.  We  will  study 
that  case  in  detail  in  Appendix  B. 
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Appendix  B:  Fast  Computation  of  the  Coefficient  of  Potential  Matrix 


The  following  summary  of  the  method  for  obtaining  the  coefficient  of  potential  matrix  is 
done  with  reference  to  [5]  and  [9], 

Recall  the  coefficient  of  potential  matrix  has  elements  py,  which  gives  the  potential  V,  on 
rectangle  i  due  to  a  charge  Q,  placed  on  rectangle  j,  is  defined  as 


Pu=—\  j  Gedvjdvt 

x;  x;  J  v;  J  V,-  J 


VVj  Jvi  Jvj 


Because  the  thickness  is  very  small  relative  to  the  surface  area,  we  ignore  the  thickness  of 
rectangles  i  and  j  and  integrate  over  the  surface,  which  gives 


Pij=  —  f  f  G  dsJdsi . 

C  C  JS;  JS;  J 


S.S  .  Js‘  JsJ 


From  Appendix  A  we  know  that  Ge  is  of  the  fonn  Ge  =  A  ( x;  x')  Y  (  v;  y'  )  Z'  (z;  z' ) .  Plugging  in 

( m7V  ^  Tr  ( rm  \  ( m7t  f  rm  .  _ 

X  =  cos  - x  ,  7  =  cos  — y  I  and  Z  =Z(z;z)cos  - x  cos  — y  ,  we  can  re-write  Ge 


as 

( mn  }\  ( nn  ^ 

cos  - x  cos  —  y  Zmn(z;z)\.  (B.3) 

\  a  J  V  b  J  J 

where  Zmn(z;z')  is  found  from  equation  (A.30).  The  set  of  rectangles  i  and  j  and  their 

coordinates  is  given  in  Figure  3 1 . 


G.-ZZ 


X  cos 
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Figure  31  -  Coordinates  for  rectangles  i  and  j  (based  on  [5]) 


Performing  the  integration  from  equation  Figure  31  on  the  set  of  contacts  for  rectangles  i  and  j 
shown  in  Figure  31,  ignoring  the  m  =  n  =  0  case,  gives 


p.-YL 


sin  nut—  -sin  nut—  sin  nut  -sin  nut 


sin  wt  —  -  sin  rut  —  sin  rut  —  -  sin  rut  — 


hi  b\ 


b,-b, 


The  sinusoidal  terms  in  equation  (B.4)  are  then  multiplied  out  and  we  use  the  identity  shown  in 
equation  (B.5)  to  put  equation  (B.4)  into  the  fonn  of  a  discrete  cosine  transform  (DCT). 


at  L1 


sin  mTt  — -  sin  nut—  =—  cos  nut 


a:  -  a , 

1  J 


-cos  mTt 


a,  +  a . 

1  J 


As  in  [5]  and  [9],  this  allows  us  to  re-cast  equation  (B.4)  into  a  sum  of  64  terms  of  the  form 


114 


cos  nut 


an  —  a3,4 


cos  nn 


.  ^1,2  —  ^3,4 


This  is  of  the  same  form  as  a  2-D  DCT.  By  inspection  of  equation  (B.6),  it  is  apparent  that  for 
the  case  m  =  0,  n  0  or  m  0,  n  =  0 ,  the  fonn  of  equation  (B.6)  becomes  the  form  of  a  1-D  DCT. 
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Equation  (B.6)  is  not  solvable  due  to  the  infinite  upper  bound  on  the  DCT.  To  make 
equation  (B.6)  solvable  with  a  DCT,  we  need  to  discretize  the  substrate  in  the  x  and  y  directions 
and  represent  the  coordinates  of  the  contacts  as  ratios  of  integers  as  in  [5]  and  [9]. 


=  P^\=(h_ 
"  0 


a  P  b 

This  allows  equation  (B.6)  to  be  re-written  as 


(B.7) 


p- 1  2-1 

YYz 

mn 

m= 0 n= 0 


r 


cos 


mn  - 


P\,2  ±  Px- 


cos 


vin  - 


<h,2±<h< 


(B.8) 


v  *  y  v  0 

The  variables  a  and  b  in  equation  (B.7)  are  the  x  and  y  dimensions  of  the  substrate.  They 
do  not  need  to  be  the  exact  value,  just  large  enough  for  convergence.  The  inductors  we  studied 
are  on  the  order  of  a  few  hundred  microns  in  diameter  and  a =b =1024  um  is  adequate  for 
convergence.  The  values  of  P  and  Q  determine  the  number  of  iterations  over  which  the  infinite 
summations  are  made,  or  in  other  words,  the  size  of  the  DCT. 

Equation  (B.8)  is  not  exactly  a  DCT  because  of  the  ±  terms  in  the  numerator  which  have 
an  upper  limit  up  to  2 P  and  2 Q.  To  truly  be  a  DCT,  the  indices  in  the  numerator  of  the  cosine 
arguments  in  the  DCT  must  not  be  greater  than  the  upper  limits,  P- 1  and  Q- 1  in  this  case.  The 
following  properties  of  symmetry  of  the  DCT  allow  us  to  get  around  this.  We  define  Zpq  as 


p- 1  e-1 

z  =  y  Vz 

pq  mn 

m= 0  n= 0 


7  =7  =7  =7 

^ P,q  ^ 2P-p,q  ^ p,2Q~q  ^ 2P-p,2Q-q 


cos 

p  } 

mn  — 

cos 

f  q) 

nn  — 

l  PJ 

l 

(B.9) 
(B.  10) 


Equation  (B.9)  is  in  the  exact  form  of  a  P  by  Q  DCT.  Thus,  if  the  terms  pX2±piA  and  qx  2  ±  q3  4 

are  greater  than  P- 1  or  Q- 1,  respectively,  we  can  flip  the  values  across  P  and  Q  using  the 
identities  in  equation  (B.10). 

For  the  case  m  =  n  =  0,  the  portion  of  Ge  that  is  a  function  of  x  and  v  drops  out  because  it 
is  equal  to  1  for  m  =  n  =  0.  Thus,  equation  (B.2)  becomes 

p7=  —  f  f  Zoo dsjdsi.  (B.  11) 

SiS.Js‘JsJ 

If  we  perfonn  this  integration,  we  find  that  p  *'  is  simply  the  Green’s  function  for  in  =  n  =  0 


since  the  integration  over  the  contacts  is  in  x  and  y  while  Zoo  is  a  function  of  z  only. 


00  /^*00  /  .  t\  r7  /  .  f\ 

Pi,  =Ge  (z;z)  =  Z00(z;z  ) 


(B.  12) 
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To  find  Zoo,  we  inspect  equation  (A.  18)  with  tmn  =  0.  The  solution  to  Z0o  is  of  the  form 


z;0  =  k(c;z+d;), 

where  the  superscript  u  indicates  z  >z  as  in  Appendix  A.  For  z'  <  z ,  we  have 

z'0=^(c;z'+jd;). 


(B.  13) 


(B- 14) 


As  in  [10],  we  find  K  by  integrating  equation  (A.  18)  for  z  =  z'  and  applying  the  boundary 
condition  that  the  potential  be  continuous.  Doing  this  gives 


K  = 


1 


absC. 


u,l 


(B.  15) 


We  choose  the  superscript  u  or  /  in  equation  (B.  15)  depending  upon  the  location  of  z  with  respect 
to  z'  as  done  previously. 

To  find  the  values  for  C"’z  and  Dl‘J  we  use  a  matrix  [A]  as  in  equation  (A.9),  but  [A]  is 


defined  differently  for  the  case  m  =  n  =  0.  If  we  apply  the  boundary  conditions  of  continuous 
potential  and  continuous  hxD  at  the  interface,  as  well  as  Z00(0;z')  =  0  due  to  the  ground  plane 


and  77  •  E  =  0  at  the  top  of  the  substrate,  we  find  [A]  as 


M- 


'k+l 


f  e  ^ 
1— ^ 


d„  1 


'k+ 1  J 


(B.  16) 


We  are  free  to  set  C;  =  1  and  we  must  set  D;  =  0  because  of  the  ground  plane  at  z  =  0.  Thus 
far,  we  have  defined  py  for  all  values  of  in  and  n.  To  summarize  the  procedure  from  [5]  and  [9], 
the  steps  for  fast  computation  of  [f]  are  given  below. 


1 .  Input  the  substrate  data.  This  includes  the  thickness  and  dielectric  constants  of 
the  dielectric  layers  of  the  substrate. 

2.  Compute  Zmn  as  outlined  in  Appendix  A  and  Z0 o  as  outlined  in  Appendix  B. 

3.  Compute  the  DCT  of  Zmn  for  all  cases  of  m  and  n  except  for  m=n=0. 

4.  Read  in  the  x  and  y  coordinates  of  the  comers  of  the  rectangles  comprising  the 
current  geometry  and  compute  the  equivalent  integer  representation  as  given 
by  equation  (B.7). 
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5.  Compute  equation  (B.4)  using  equation  (B.8).  Recall  that  for  m=n=0,  use 
Pij  =Ge  (z,z)  =  Z00(z;z  ) . 

6.  Generate  [P]  by  computing  all  possible  combinations  of  rectangles.  For 
example,  if  we  have  an  inductor  of  20  rectangles,  [P]  will  be  20  x  20. 

7.  Steps  1-3  only  need  to  be  done  for  a  given  substrate.  If  [ P ]  needs  to  be 
recomputed  for  a  different  geometry,  go  to  step  4. 
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Appendix  C:  Mutual  Inductance  Calculation 

The  magnetostatic  Green’s  function,  Gm,  is  used  to  calculate  the  partial  inductance  matrix 
\M\.  We  make  the  assumption  that  the  dielectric  substrate  under  consideration  has  low  loss, 
which  leads  to,  from  equation  (2.3) 

V2A  =  ~/uJ .  (C.l) 

By  inverting  equation  (C.l)  and  setting  J  to  be  a  dirac  delta  function,  we  get  Gm.  Since  //  =  //„ 
for  all  layers,  Gm  is 


Mo 

AnR 


As  discussed  in  Chapter  2,  we  find  My  from  Gm  as 


My  =  - — - 

4  nap. 


1111, ITT Tda‘dar 


j  j  \r  -  r 


(C.2) 


(C.3) 


In  this  appendix,  we  will  review  the  techniques  of  [2], [6],  and  [7]  to  solve  equation  (C.3). 

We  begin  with  the  diagonal  tenns  of  \M\,  also  known  as  the  self  inductance  terms.  For  a 
rectangle  of  thickness  t,  width  W,  and  length  L,  we  get  Mu  to  be 


M,  () nH )  =  2  x  KT4  x  L(jum)  x  fin  f - 2L{MIU) - 1  +  0.5I  (C.4) 

The  calculation  of  the  off-diagonal  terms  of  [M],  also  called  the  mutual  inductance  terms,  is 
more  complicated.  The  calculation  of  the  mutual  inductance  has  two  parts:  the  4  dimensional 
integral  and  the  dot  product. 

The  dot  product  is  quite  simple.  If  the  rectangles  are  oriented  in  such  a  way  that  the 
lengths  are  perpendicular,  then  t//;  •  dlj  is  zero.  If  the  lengths  are  in  parallel  or  in  series,  then 


•  dlj  =  1  for  current  in  the  same  direction.  The  dot  product  dlt  •  dlj  =  -1  if  the  current  is  in  the 
opposite  direction. 

The  calculation  of  the  4  dimensional  integral  is  done  in  one  of  two  ways,  depending  upon 
the  relation  of  the  center-to-center  distance  D  and  length  of  the  rectangles  L.  First,  we  treat  the 
rectangles  as  filaments  with  no  width  and  the  same  length  L  and  separated  by  a  distance  d.  The 
solution  to  the  integral  of  equation  (C.3)  for  this  case  is 
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The  criteria  for  deciding  when  to  use  the  til  amen  tal  approximation  of  equation  (C.5)  or  the 
GMD/AMD  approximation  of  equation  (C.9)  is  based  upon  the  size  of  d  in  relation  to  the  length 
of  the  rectangles.  If  d  <  .5 L ,  the  widths  of  the  rectangles  cannot  be  ignored  and  we  use  the 
GMD/AMD  approximation.  If  d  >  .5  L  ,  we  use  the  til  amen  tal  approximation. 
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In  general,  the  two  rectangles  will  not  be  the  same  size  and  will  sometimes  have  an 
offset.  [2]  gives  a  method  for  calculating  the  mutual  inductance  of  two  rectangles  that  are  offset 
and  not  equal  in  length.  There  are  four  different  possible  orientations  for  two  rectangles  of 
arbitrary  length  and  offset.  The  first  is  given  in  Figure  33. 

I 


d 


Figure  32  -  Unequal  parallel  filaments  with  offset 


The  mutual  inductance  M'  of  the  rectangles  with  this  configuration  is  given  as 

M'  =  \{Ml+m+s+Ms)~(Ml+s+Mm+s),  (C.10) 

where  the  subscripts  are  the  values  for  L  to  be  used  in  equation  (C.5)  and  (C.9).  Equation  (C.10) 
is  for  the  case  when  the  offset  is  positive,  i.e.  the  filaments  do  not  overlap.  If  the  filaments  do 
overlap,  then  we  have 

M'  =  +M  (C.1 1) 

If  d  and  5  go  to  zero  in  Figure  32,  we  get  two  rectangles  in  series.  For  the  case  when  d  and  8  go 
to  zero  we  use 

(C.  12) 

A  variation  of  Figure  32  is  when  there  is  a  full  overlap  of  the  rectangles  as  shown  in  Figure  33. 
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Figure  33  -  Unequal  parallel  filaments  with  full  overlap 
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The  mutual  inductance  for  the  case  of  Figure  33  is 


(C.  13) 


Lastly,  there  is  the  case  when  both  filaments  are  unequal  in  length  but  have  no  offset  as  shown  in 
Figure  34. 


d 

m 

Figure  34  -  Unequal  parallel  filaments  with  no  offset 

The  mutual  inductance  for  this  case  is  found  as 

(C.14) 
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